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Abstract 

The stability of some spatial asymmetric games is discussed. Both 
linear and nonlinear asymptotic stability of asymmetric hawk-dove 
and prisoner's dilemma are studied. Telegraph reaction diffusion equa- 
tions for the asymmetric spatial games are presented. Asymmetric 
game of parental investment is studied in the presence of both ordi- 
nary and cross diffusions. 

1 Introduction 

In asymmetric games [Hofbauer and Sigmund 1998], different players have 
different strategies and different payoffs. In reality, most games are asym- 
metric e.g. battle of the sexes and owners-intruders games. The differential 
equations of the asymmetrical games are 

§ = » [(Bp), - qBp] W 

where A(B) is the payoff matrix of the first (second) player, and Pi(qi) is 
the fraction of adopters of the strategy i in the first (second) population, 
respectively. 
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The problem of Turing instability (diffusion induced instability) [Okubo 
f 980] for symmetric games have been already studied [Cressman and Vickers 
1997]. In this case, the standard replicator equation for the symmetric game 
is given by 

^=p i [(Ap) t ~ P Ap] + DV 2 p l (2) 

It is known [Okubo 1980] that including spatial effects may significantly 
change the stability of equilibrium points. Also spatial effects are crucial 
in many biological phenomena. Some authors [Boerlijst and Hogeweg 1995] 
have studied general spatial hypercycle systems. Therefore studying spatial 
asymmetric games is an important problem. Due to the difficulty in defin- 
ing evolutionarily stable strategy (ESS) in asymmetric games [Hofbauer and 
Sigmund 1998], we will concentrate on asymptotically stable strategies. The 
equations of spatial asymmetric games are [Hofbauer et al 1997] 

$= Pi [(Aq) i -pAq]+D 1 V 2 p i , 

%= gi [(Bp)i - qBp] + D 2 V\. { ' 

In this paper, we will attempt to answer the following questions: 

1. Given an asymptotically stable solution to the system (1), does Turing 
instability exist for the corresponding spatial game (3)? 

2. Given an asymptotically unstable solution to the system (1), can dif- 
fusion stabilize it? 

3. Given an asymptotically linearly stable solution to the system (1), is it 
nonlinearly stable? 

Our typical examples will be the asymmetric hawk-dove (AHD) and the 
asymmetric prisoner's dilemma (APD) games. 

The paper is organized as follows: In section 2, the asymmetric hawk- 
dove game is studied. Conditions for Turing stability and nonlinear finite 
amplitude instability are derived. In section 3, The asymmetric prisoner's 
dilemma is presented. Telegraph reaction diffusion equation is applied for 
stability analysis of the asymmetric prisoner's dilemma game in section 4. In 
section 5, an asymmetric game of parental investment will be studied. Some 
conclusions are summarized in section 6. 
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2 Asymmetric hawk-dove game 



In this case the possible strategies are hawk (H) and dove (D), and the payoff 
matrices A and B in Eq. (f) become 



|(fi-ci) Vi 
a 



B 



2 {v 2 - c 2 ; 




V2 
2 



(4) 



where q > fj, i = 1,2. The corresponding partial differential equations for 
the spatial AHD are 



f t =D 1 V 2 p + ±p(l-p)(v 1 -c 1 q), 
f t =D 2 V 2 q + \q(l-q)(v 2 -c 2 p), 



(5) 



where p (q) is the fraction of hawks in the population of the first (second) 
player. It is direct to see that the solution p — 1, q = is linearly asymptot- 
ically stable solution for the system (5), without diffusion (D 1 — D 2 — 0). 

The first question is about Turing (diffusion induced) instability [Okubo 
1980]. It occurs if the following system 



fli& + /i(p,g), 



% = D 2 H + f 2 (p,q), 



dt 1 dx 

dg n 9 

dt ~ 1J2 dx 2 

has an equilibrium solution (p ss , q ss ) which is stable if Di 
corresponding linearized system 



p 



satisfies the condition 



p ss + e(x,t), q 



q ss + r](x,t), 
+ due + a i2 rj, 



at ^ 1 d% 

f = D 2 y + a 21 e + a 22 r], 



(6) 

D 2 = 0, and the 



(7) 



H(k 2 ) = DiD 2 k A - {DXCL22 + D 2 a u )k 2 + a ±1 a 22 - a 12 a 21 < (8) 

In this case diffusion will destabilize the solution (p ss ,q ss ). This is Turing 
instability. 

Applying the above procedure to the spatial AHD one gets: 



Proposition (1): The equilibrium solution p 
game is Turing stable. 



1, q = of the AHD 
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The second question to be discussed is: Can diffusion stabilize an unstable 
solution of the AHD game? Consider the internal solution p = V2/C2, q = 
v\jc\. It is asymptotically unstable solution to the diffusionless case [D\ = 
D 2 = 0). Including diffusion and linearizing around p = v 2 /c 2 , q = V\/c\, 
and assuming the following boundary conditions: 

p=f 2 +e(x,t), q=f i +r ] (x,t), 
e(0,i)=e(l,f) = 0, 77(0, *) = 77(1, t) = 0, 1 ) 

then 



Proposition (2): The interior solution p = v 2 /c 2 , q = Vi/ci, with the 
boundary conditions (9) is linearly asymptotically stable if 



Dl D 2 «'-^(l-*)(l-*)>0. (10) 
4 ci c 2 



Linear stability analysis is useful if the perturbations of equilibrium are 
infinitesimally small. This is not always the case in biological systems. In 
this case one has to study finite amplitude instability (FAI) [Stuart 1989] of 
the equilibrium solution. In the following, we generalize the work of Stuart 
to the two species case. Consider the following equation 

Linearizing around the solution = 0, i.e. let 



9(x, t) = v(x, t), 

linearize in v, then 

t = + W(O), Letv = <P(x) expert) h9) 

a<j> = r + r (0)0, 0(0) = 0(1) = 0. [LA) 

Set = Y^i^i sm(irlx), then bifurcation points are given by /'(0) = (nl) 2 , I = 
1,2,.... Studying the stability of the first bifurcation point / = 1 using 
Matkowsky two-time nonlinear stability analysis, one defines A = a — k, 
set A = 7r 2 + Ao£ 2 , e is a small parameter. Decompose the time into fast t' 
and slow r, then 
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Expand u in powers of e 2 , then u ~ ev\ + e 3 v 3 + . . ., (notice that /"(0) = 0), 
then substituting in Eq. (12), one gets 

vi(x, t', t) — oz(t) sin(7r/a;) exp(l — l 2 )n 2 t'. 
i 

Substituting into the cubic term and setting the constant term in t' equal to 
zero, one finally gets 

^ = \ a 1 + b(a 1 )\ 
t = 4,V'(0)^. 

D J sin 7ra; ax 

Thus a nonlinear (finite amplitude) instability arises if 



A < and M0)| > (13) 

o 

In the beginning of the section the condition f"(0) = was imposed, here 
we will assume f"(0) ^ 0. Thus we consider 

nil 

— = DV 2 u + f(u), f(0) = 0, u(0, t) = u(l, t) = 0. (14) 
The solution u = is a steady state solution, so expanding near it we set 

d d d 



u=Y,e m v m {t' ) t" ) x) ) _ = _ + e _, /'(0) = Ao + eAi. 



m=l 

Substituting in Eq. (14) and equating terms 0(e), one gets 
Vl (t', t", x)=J2 a i(t") sm(7r/a;) exp[(l - Z 2 )tt 2 ]. 

Let Ao = 7r, and consider the equation 0(e 2 ), we set the secular term (inde- 
pendent of t') equal to zero, then 

dai , fb'„, nS \ 2 r/ /(! sin 3 7ra; dx , . 

- = Al0l +(-/») al, V = (15) 



Thus the conditions for FAI are 

2A 

f(0)>0, Al < and Ql (0) > 

For systems of two partial differential equations 

= AW + f(u u u 2 ), = Z) 2 V 2 « 2 + ua), 
ui(0,*)= ui(l,*) =u 2 (0,t) = u 2 (l,t) = 0, (16) 
/(0,0) = «?(0,0) = 0. 

Expanding near the steady state solutions ui = u 2 = 0, we get 



Ml = E m =l £ m Wm(t', t", X), U 2 = Em=l ^m^', t", 
9. _9_ _|_ 9 

at at' " r fc at" ' 



~~ 9uf(0,0) ~~ ^ n + £ ^ 12 ' ^ 2 ~~ ^ 21 + £ ^ 22 ' 
01 = A*ll + £A*12, 02 = A*21 + ^22- 

After some tedious calculations, we got the following conditions for FAI in 
the system (16): 

(i) (n 2 D 1 - \ 11 )(ir 2 D 2 - fji 21 ) - \ 2mi = 0. 

(ii) A12 + k\ 22 = *f + fi 22 < 0, where k = ^ D ^ 11 ■ 

(iii) ff + K f 12 + ^f 22 = f: + g 12 + fg 22 >0. 

(iv) ai(0) > /f Al2+KA22 2 v 

where 6' is defined in Eq. (15). Applying the condition (i) to the spatial 
AHD, one gets 

(n 2 D 1 +v 1 )(7c 2 D 2 + v 2 ) < 0, 
which is not possible, thus we find: 

Proposition (3): There is no finite amplitude instability for the solution 
p = 1, q = of the AHD game. 
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3 Asymmetric prisoner's dilemma game 



In the prisoner's dilemma game, the possible strategies are cooperate (C) 
and defect (D). The payoff matrices A, B in Eq. (1) are given as follows: 



" Ri 


Si ' 


, B = 


' R 2 


s 2 ' 


. T i 


Pi 




P 2 



(17) 



such that 2Ri > Tj + Si, and Tj > Ri > Pi > Si, where i — 1, 2. The 
dynamical equations for the spatial asymmetric prisoner's dilemma game 
(spatial APD) are: 

^ = D 1 W 2 u 1 + 
Wl (l - «i)[-(Pi - 50 + U2 (P! - S 1 - T x + R ± )}, 

dm = D 2 V 2 u 2 + 
u 2 {\ - u 2 )[-(P 2 - S 2 ) + Ul (P 2 -S 2 -T 2 + R 2 )}, 

where ^1(^2) is the fraction of cooperators in the first (second) players pop- 
ulation. The solution u\ = u 2 = is linearly asymptotically stable. It 
represents the always defect strategy. 

Two questions arise the first is: can diffusion stabilize the cooperation 
solution ui — u 2 — 1? And does the always defect solution have (FAI) 
nonlinear instability? Using the techniques of the previous section we get: 

Proposition (4): 

(i) If Dili 2 > Ti — Ri, then the cooperation solution is linearly asymptoti- 
cally stable. 

(ii) The always defect solution does not have FAI. 



4 Telegraph reaction diffusion in spatial asym- 
metric games 

The standard spatial games depend on the familiar reaction-diffusion equa- 
tion 

9u^ = D Sh^ + m (lg) 
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A basic weakness in this equation is that the flux j reacts simultaneously to 
the gradient of u consequently an unbounded propagation speed is allowed. 
This manifests itself in many solutions to Eq. (1) e.g. (if / = 0), then 

1 ^£ 

u(x, t) = e 4Dt ) u (x, 0) = 5(x) i.e u(x, t) > OVx, Vt > 0. 

This is unrealistic specially in biological and economical systems, where it is 
known that propagation speeds are typically small. To rectify this weakness, 
Fick's law is replaced by 

dj du 

and the resulting telegraph diffusion equation becomes 

d 2 u du „ d 2 u 

r 1 = D . 

dt 2 dt dx 2 

The corresponding telegraph reaction diffusion (TRD) is 

The time constant r can be related to the memory effect of the flux j as a 
function of the distribution u as follows: Assume that [Compte and Metzle 
1997] 



j(x,t) = -J*K(t-t')^^dt', (21) 



^du(x,t') 

jyx,i) = - I j\yi - 
Jo 

hence 

j + T % = « - H (r 9 -^ + K (t - i0 ) | df . 

This equation is equivalent to the telegraph equation if 

D ,-f N 
K(t) = -exp(—). 
r r 

This lends further support that TRD is more suitable for economic and bi- 
ological systems than the ordinary diffusion equation since e.g. it is known 
that we take our decisions according to our previous experiences, so memory 
effects are quite relevant. Further evidence comes from the work of Chopard 
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and Droz [Chopard and Droz 1991], where they have shown that starting 
from discrete time and space then the continuum limit does not give the 
standard reaction diffusion but the telegraph one. 

Since it is known that new technologies, habits etc... takes time to spread, 
we believe that TRD equation is more relevant than the ordinary diffusion 
equation in modelling economic and biological systems [Ahmed et al 2001]. 

Now we apply TRD to spatial APD game. The TRD for a system of 
equations are [Hadeler 1998] 



(22) 



' at 2 ^ at 1 at d Uj 
DV 2 Ui + fi(u u u 2 , . . .,u n ), 

hence applying it to the APD (18), we get 

T § & + ! &-T*t(l-2u 1 )[-{P 1 -S 1 ) + 

u 2 (Pi - Si - 7\ + i?i)] = £i W+ 
ui(l - ui)[-(Pi - Si) + u 2 (P 1 - S 1 - Ti + i?i)], 

^ + 3 -aT ~ - 2«2)[-(P 2 - S 2 ) + 

Ui(P 2 -S 2 -T 2 + R 2 )} = D 2 V 2 u 2 + 
u 2 (l - u 2 )[-(P 2 - S 2 ) + Ul (P 2 -S 2 -T 2 + R 2 )]. 

The following question arises: Can diffusion stabilize the cooperation so- 
lution ui = u 2 = 1? Using the techniques of the second section, we get 



(23) 



Proposition (5): If the following conditions are satisfied 

D i n 2 >T i -R i , 1 > riTi-Ri), 
4r(Avr 2 -Ti + Ri)<[-1 + r(T, - R,)} 2 , i = 1,2, 

then the cooperation solution is linearly asymptotically stable. 

Proof. Assume that 

u± — 1 — E\ exp(o"it) sin(7rx), u 2 — 1 — e 2 exp(a 2 t) sm(nx). 
Substituting one gets 

*i = ±[(-1+t(T 1 -R 1 ))± 

+ r(Ti - i?i)) 2 - 4r (R, - T x + D 1 tt 2 )}, 
a 2 = l T [{-l+T(T 2 -R 2 ))± 

+ r(T 2 - R 2 )f - At (R 2 - T 2 + D 2 ^)}. 



(24) 
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Stability requires that the real part of a i: % — 1, 2 is negative. The first two 
conditions in the proposition guarantee this requirement. Furthermore since 
Ui, % — 1,2 are real and nonnegative, then the term under the square root 
should be nonnegative. The third condition of the proposition guarantees 
Ui > 0. 

Thus the conditions for cooperation stability for TRD are more stringent 
than those for ordinary diffusion (c.f. proposition (4)). 

5 Asymmetric game of parental investment 

Parents are faced with the decision whether to care for the offsprings or to 
desert. A model has been given for this asymmetric game [Krebs and Davies 
2000]. Let po, pi, p 2 , be the probabilities of survival of offsprings which are 
not cared for, cared for by a single parent and cared for by both parents, 
respectively, then p < p\ < p 2 . A deserting male has a chance q of mating 
again while a caring (deserting) female has Wi(w 2 ) offsprings. The payoff 
matrices for male (female) corresponding to the strategies C (care) or D 
(desert) are denoted by A(B), and given by 

[ iuipi(l + q) w 2 po{l + q) 

The spatial asymmetric equations for the above game are: 

^ = u[wip 2 v + tu 2 Pi(l -v) - wip 2 uv - w 2 p 1 u{\ — v) — 
w lPl (l + q)v(l -u)- w 2 p (l + q)(l- u)(l - v)] + + D 12 ^, 

|| = v[wip 2 u + wip 2 (l - u) - W\p 2 uv - w 2 piu(l - v) — 
w lPl v(l u) w m {\ -u){\- v)\ + D 2 + £> 21 0. 

In this system, we introduced both ordinary and cross diffusion. Cross dif- 
fusion is the diffusion of one type of species due to the presence of another 
[Okubo 1980]. This phenomena is abundant in nature e.g. predator-prey sys- 
tems where the predator diffuses towards the regions where the prey is more 
abundant. On the other hand the prey tries to avoid predators by diffusing 
away from it. Another area of application is in epidemics where susceptible 
individuals try to avoid infected ones. 

Here we will see that ordinary diffusion is unable of destabilizing the 
diffusionless ESS: 



B 



wip 2 Wipi 
w 2 pi w 2 p 



(25) 
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(i) ESS1 where both male and female desert i.e. (u — 0, v — 0). It requires 
w 2 p > W\P\ and p (l + q) > p 1 . 

(ii) ESS2 where male cares and female desert i.e. (u — 1, v — 0). It requires 
w 2 pi > and p (l + q) < p x . 

(iii) ESS3 where female cares and male desert i.e. (u — 0, v — 1). It requires 
WiPi > w 2 p and + g) > p 2 - 



[iv) ESS4 where both male and female care i.e. (u 
Wip 2 > w 2 pi and + q) < p 2 . 



1, v — 1). It requires 



Following steps similar to the previous games, we get 

Proposition (6): The solution ESSi, % = 1,2,3,4 is destabilized if the 
following condition is satisfied: 



where 



D 12 D 21 ir 4 > (D l7 i 2 - ai ){D 2 Tx 2 - k) 



ai = w 2 pi - w 2 po(l + q), h = wipi - w 2 po, 

a 2 = -u! 2 pi + u> 2 R)(l + g), h = Wip 2 ~ w 2 p u 

a 3 = w x p 2 - + q), bi = -Wipi + w 2 p , 

a 4 = -wtfz + wjp^i + q)^ ^ = w 2 pi - w x p 2 . 



(26) 



Notice that all a«, 6j, % = 1,2,3,4 are negative hence ordinary diffusion 
cannot destabilize the ESS. 

Applying the above procedure to the battle of the sexes [Schuster and 
Sigmund 1981] where the female has two strategies coy or fast while the 
male can be either faithful or philanderer. The male (female) payoff matrix 
is A(B) 







-10 



-2 



B = 





3 



hence the spatial battle of the sexes equations are 



du 
dt 



at 



u(l-u)(-10 + 12v) + D i y 
l-u)(5-8«) + D 2 0. 



(27) 
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Proposition (7): Diffusion stabilizes the internal equilibrium of the sys- 
tem (27). 



Proof. There is a unique internal homogeneous equilibrium solution E = 
(5/8,5/6), which (for the diffusionless case) is stable but not asymptotically 
stable. Substituting with 

u — — h e exp(crt) sin ttx, v = — h q exp(crt) sin nx, 1 > x > 0, 
8 6 

in Eq. (27), and linearizing in e, (, one gets 

(a + D^)e = ^C, {* + D^ 2 )C = ~e. 

Hence 

25 

[( ( 7 + J D 1 7T 2 )( ( 7 + J D 2 7r 2 ) + y ]C = 0, 

25 

(a + D 1 n 2 )(a + D 2 n 2 ) + — = 0, 

o 

95 

a 2 + ^(DiTr 2 + D 2 n 2 ) + (4 + D i D 2^) = 0. 

8 

Therefore the real part of a is negative, then the internal equilibrium (in the 
presence of diffusion) is asymptotically stable. It is clear that if the diffusion 
coefficients are set equal to zero (Di — D 2 — 0), then one regains the stability 
but not asymptotic stability. This completes the proof. 



i.e. 



then 



6 Conclusions 

Based on replicator equations, a mathematical approach for the analysis of 
spatial asymmetric games is introduced. Some questions regarding spatial 
stability for asymmetric hawk-dove and the asymmetric prisoner's dilemma 
(APD) games are answered. Telegraph reaction diffusion equation is applied 
for stability analysis of the asymmetric prisoner's dilemma game. Asymmet- 
ric game of parental investment is studied in the presence of both ordinary 
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and cross diffusions. Ordinary diffusion cannot destabilize the ESS for this 
game. Conditions for destabilizing the ESS are given in the case of cross 
diffusion. 
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